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ABSTRACT 


Optical  scintillation  on  folded  paths  has  been  investi¬ 
gated  for  the  case  of  spherical  waves  and  a  plane  mirror  as 
a  folding  target.  A  criterion  for  the  statistical  indepen¬ 
dence  of  the  two  parts  of  the  folded  path  was  formulated, 
utilizing  the  path  integral  technique  and  the  Von  Karman 
spectrum.  Employing  the  same  approach,  an  expression 
for  the  scintillation  strength  (the  log-intensity  vari¬ 
ance)  on  folded  paths  has  been  developed.  It  has  been 
shown  that  the  scintillation  strength  on  folded  paths 
earn  be  expressed  as  a  product  of  the  one-way  scintilla¬ 
tion  strength  and  a  function  of  the  distance  source-detector 
in  units  of  the  outer  scale.  The  analytical  prediction  3hows 
that  the  ratio  of  the  log  intensity  variances  on  the  folded 
paths  to  the  one  on  the  one-way  path,  goes  from  ,27^T^  for 


exact  folding  (the  distance  source-detector ip  zero) ,  to 
-^^^when  the  detector  is  far  away  from  the  source.  The 
theoretical  predictions  have  been  corroborated  by  preliminary 


field  experiments. 
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I.  INTRODUCTION 


Intensity  fluctuations  of  visible  or  near-visible  sources, 
while  viewed  through  the  atmosphere  (scintillation) ,  play  a 
major  role  in  a  large  variety  of  situations  such  as  star 
observations,  optical  communication  channels,  Lidars,  power 
transmission,  etc.  It  has  been  understood  for  a  very  long 
time  that  the  stars'  twinkling  is  caused  by  the  inhomogeneity 
of  the  atmospheric  index  of  refraction.  It  is  only  in  the 
last  thirty  years,  and  especially  since  the  first  lasers  were 
built  about  twenty  years  ago,  that  this  intuitive  understand¬ 
ing  has  been  put  into  a  rigorous  quantitative  model.  Some 
of  the  better  known  contributors  to  this  progress  are 
Tatarski  [1,2],  Fried  [3],  Lutomirski  fund  Yura  [4],  Ishimaru 
[5],  Lee  and  Harp  [6],  Clifford  [7],  Fante  [8,9]  and  Dashen 
[10] .  As  a  result  of  their  work  and  many  others  it  is  now 
believed  that  in  most  of  the  situations  the  processes  involved 
are  fairly  well  understood  and  the  theoretical  predictions 
are  reasonably  close  to  experimental  findings. 

Among  the  problems  which  are  less  understood  and  for  which 
existing  theory  either  does  not  give  reasonable  predictions 
or  has  not  been  confirmed  experimentally,  is  the  problem  of 
scintillation  on  folded  paths.  In  general,  this  descriptive 
name  stands  for  a  class  of  situations  in  which  EM  fields 
propagate  in  the  atmosphere  from  a  source  to  a  target  and 
beck  to  a  receiver  (detector)  in  the  vicinity  of  or  at  the 
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same  location  of  the  source.  Typical  situations  which  belong 
to  this  class  are  optical  radars,  communication  links  (ground- 
satellite/mirror-ground)  ,  remote  monitoring  of  atmospheric 
pollutants,  and  others. 

The  class  of  point-ahead  problems  such  as  pointer-tracker 
(with  fast  moving  targets) ,  power  transmission  from  space  to 
earth  (which  follow  a  beacon  signal  from  earth  to  space) 
shares  a  major  question  with  the  previous  class,  that  is  to 
say,  the  correlation  between  the  two  parts  of  the  path. 

In  this  paper  we  are  concerned  with  scintillation  caused 
by  propagation  along  a  folded  path  when  the  target  is  a 
specular  reflector  (either  a  plane  mirror  or  retroref lector) . 
Transition  cases  in  which  the  two  parts  of  the  path  are  par¬ 
tially  correlated  are  included.  To  the  best  of  our  knowledge, 
only  a  small  number  of  works  in  the  open  literature  address 
this  problem.  We  will  compare  the  results  of  these  works 
with  ours  at  the  related  parts  of  our  paper. 

In  the  second  part  of  this  paper  (discussion)  we  describe 
some  accepted  formulas,  physical  pictures  and  approximations, 
which  describe  the  atmospheric  turbulence  and  the  propagation 
of  EM  fields  in  it.  In  the  third  part  we  develop  an  expression 
for  the  correlation  between  the  two  parts  of  the  path.  In 
the  fourth  we  give  the  expressions  for  the  scintillation 
strength  (the  normalized  variance  of  the  log-intensity)  for 
a  plane  mirror  and  retroreflector  targets.  The  fifth  part 
describes  the  experimental  part  of  our  work  (system  descrip¬ 
tion  and  results),  and  a  summary  constitutes  the  sixth  part. 
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II.  BACKGROUND 


A.  THE  SCALAR  EQUATION  OF  PROPAGATION 

The  effect  of  the  atmosphere  on  the  EM  fields  propagation 
manifests  itself  through  the  Maxwell  equations.  The  atmos¬ 
phere  is  considered  to  have  no  free  charges  (hence  zero  con¬ 
ductivity  and  unit  magnetic  permeability) .  It  is  also 
assumed,  based  on  experimental  findings  [8],  that  the  time 
scale  of  changes  of  the  atmospheric  index  of  refraction  is 
much  longer  than  the  time  scale  of  the  EM  field  fluctuations 
in  the  visible  and  near  visible  band.  This  assumption  allows 
us  to  suppress  the  time  dependent  part  of  the  EM  field  in 
Maxwell's  equations.  Based  on  the  above  conditions,  we  find 
that  the  vector  equation  for  the  EM  field  propagation  is 

V2E(r)  +  k2n2(r)E(r)  ?  V  (V-E(r)  )  (II-l) 


where : 


k  =  a)/c  is  the  EM  field  wave  number 
n2(r)  -  e  (r) 

e  *  the  dielectric  constant 
n  »  index  of  refraction. 

It  was  found  analytically  [11,12]  and  confirmed  in  experi¬ 
ments  [13] ,  that  the  depolarization  of  EM  fields  in  the 
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forward  direction  due  to  atmospheric  turbulence  (variations 
in  n(r))  is  negligible.  Thus  the  term,  on  the  right  hand 
side  of  Eq.  (II-l)  which  represents  the  depolarization 
(through  mixed  derivatives)  can  be  neglected  [8],  and  we  get 
the  scalar  equation  of  propagation 

V2E (r)  +  k2n2(r)E(r)  =  0  (II-2) 

B .  ATMOSPHERIC  TURBULENCE 

As  we  can  see  from  Eq.  (II-2) ,  the  effect  of  the  atmos¬ 
phere  on  the  EM  field  propagation  is  caused  by  the  fluctua¬ 
tions  of  the  index  of  refraction.  The  index  of  refraction  is 
a  very  complicated  function  of  the  atmosphere's  constituents 
(aerosols,  water  vapour  ,  etc.),  the  absorption  coefficients 
[14]  (continuous  and  discrete) ,  the  temperature  and  the  pres¬ 
sure.  Assuming  the  absorption  is  negligible  (the  wavelength 
does  not  match  any  absorption  line  or  band) ,  the  humidity 
gradient  is  negligible,  and  the  aerosols'  contribution  to 
forward  scattering  is  small,  then  the  index  of  refraction 
turns  out  to  be  a  relatively  simple  real  function  of  the 
temperature  and  the  pressure  [15] 

n(r)  *  1.0  +  Jyjy  P(r)  [1  +  -^^5]xl<T6  (II  -3) 

"—  A 

where: 

P  «  total  pressure  in  millibars 


T  -  temperature  in  degrees  Kelvin 
X  =  wavelength  in  micrometers 

Even  though  the  expression  for  the  index  of  refraction  is 
a  simple  one,  it  is  impossible  to  get  an  exact  solution  for 
the  scalar  equation  because  we  cannot  specify  the  values  of 
the  temperature  and  pressure  along  the  path.  The  only  way 
we  can  solve  this  equation  is  through  the  use  of  their  sta¬ 
tistical  properties. 

This  approach  leads  us  to  express  the  index  of  refraction 
in  the  following  form 

n (r)  *  nQ  +  n1 (r )  (11-4) 

where  nQ  is  the  average  value  of  n(r)  along  the  path  and 
n^(r)  is  the  spatial  fluctuation,  n^  is  found  by  replacing 
P(r),  T(r)  in  Eq.  (II-4)  by  PQ,  Tg  respectively,  and  it 
obviously  cannot  cause  scintillation.  In  almost  all  works 
about  the  effect  of  the  atmospheric  turbulence  nQ  is  taken 
to  be  1. 

Since  the  pressure  fluctuations  (in  the  microscale)  are 
relatively  small  116],  the  fluctuations  of  n^.  (r)  are  governed 
by  temperature  fluctuations. 

Due  to  their  small  order  of  magnitude  (less  than  one 
degree) ,  the  main  processes  through  which  temperature  inhomo- 
geneities  in  the  atmosphere  dissipate,  are  convection  on  a 
large  scale  and  molecular  dissipation  on  a  very  short  scale. 
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Therefore,  temperature  inhomogeneities  (fluctuations  in 
space)  which  have  scale  size  larger  than  a  certain  critical 
value,  are  closely  related  to  the  velocity  field  within  the 
atmosphere . 

Using  this  relationship  between  these  fields  and  the 
Kolmogorov  model  of  the  velocity  field,  Tatarski  [1]  gave  a 
statistical  description  of  the  index  of  refraction  fluctua¬ 
tions  in  the  atmosphere. 

The  Kolmogorov  [17,18]  model  is  based  on  the  following 
assumptions : 

1.  The  fluid  is  incompressible. 

2.  The  fluid  is  in  a  steady  state. 

3.  The  fluid  is  getting  a  constant  power  per-unit  mass 
from  an  outside  source. 

4.  The  Reynolds  number  of  the  flow  (onset  by  the  outside 
source)  is  much  greater  than  the  critical  Reynolds 
number  of  the  air. 

5.  Kinetic  energy  dissipation  into  heat  can  be  done 
only  on  a  very  short  length  scale. 

Since  the  Reynolds  number  of  the  original  laminar  flow 
is  larger  than  the  critical  value  for  air,  once  the  flow  is 
disturbed,  it  will  turn  into  a  turbulent  flow  with  vortices 
whose  length  scale  is  limited  only  by  the  dimensions  of  the 
flow.  This  scale  is  called  the  outer  scale  of  the  turbulence 
(Lq) ,  and  in  the  surface  layer  it  is  usually  taken  to  be  the 
height  of  the  point  above  the  ground. 
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If  the  Reynolds  number  of  the  flow  within  this  vortex  is 
larger  than  the  critical  Reynolds  number,  the  vortex  will 
break  down  into  smaller  vortices.  Since  the  Reynolds  number 
is  proportional  to  the  scale  size  of  the  flow,  the  formation 
of  smaller  vortices  reduces  the  local  Reynolds  number  within 
each  vortex  until  it  is  smaller  them,  the  critical  value,  and 
the  flow  is  stabilized.  The  smallest  vortex  size  is  called 
the  "inner  scale  of  the  turbulence"  (lQ) . 

In  the  literature  about  EM  fields  in  the  atmosphere, 
these  vortices  are  sometimes  called  eddies. 

C.  ENSEMBLE  AVERAGING,  STRUCTURE  FUNCTION  AND  TURBULENCE 

SPECTRUM 

The  averaging  process  through  which  we  get  the  moments 
of  the  EM  field  should  be  done  over  many  points.  Each  solu¬ 
tion  of  the  scalar  equation  represents  a  different  n. (r) 

•  X  — " 

field  (the  ensemble  average).  Since  the  n^(r)  field  has  a 
relatively  short  correlation  length,  averaging  over  the 
distribution  of  n^(r)  along  a  propagation  path  can  replace 
the  ensemble  average.  In  the  literature  the  ensemble  aver¬ 
age  is  denoted  by  <  >.  He  will  use  the  same  notation. 

Assuming  the  n^r)  field  is  isotropic  and  stationary, 
we  can  define  the  following  functions. 

1.  The  autocovariance  (sometimes  called  correlation 
or  covariance)  function 

2*  *  <n1(r1)n1(r2)  >  (II-5) 
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2.  The  structure  function 


*  <lni(£l)“nl(-2)1  > 


(XI-6) 


3.  The  spectrum  function  vK> 


-H*.  iK(r,-r7) 

B  ( r i » r7 )  *  ///♦n(»e  d^ 

n  “A 


Since  we  assumed  isotropic  turbulence,  the  last 
equation  can  be  written  as 


f  .  >_>  Sin  (Kp)  «f2  j». 
Br(p)  *  4tt  /  <|>n(K)  kp  K  * 


(11-7) 


in  the  last  four  equations: 


»  spatial  vectors 


ISl  ’  —2 


spatial  frequency  vector  (of  the  n^r) 
field) 


Using  the  Kolmogorov  model,  the  spectrum  function  is 


given  by 


VK) 


for  2tt/Lq  <  K  <  2ir/£0  (II-8) 
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2 

This  interval  is  called  the  inertial  range.  C 

n 

represents  the  strength  of  the  turbulence  and  it  is 
proportional  to  the  variance  of  n^  at  a  point  ri . 

In  most  early  works,  this  spectrum  was  used  (by  way  of 
approximation)  in  the  range  0  <_  K  £  This  procedure  as 
pointed  out  by  Yura  [19] ,  can  cause  serious  deviation  in  the 
outcome  of  some  calculations.  Thus  a  somewhat  different 
spectrum  is  now  being  used,  namely  the  modified  Von  Karman 
spectrum. 

2  -<KV2 

.033C„  e  u 

4>  <K)  -  - S -  (II-9) 

<K2  +  L“2)11/6 

Hill  and  Clifford  [20]  showed  that  a  high  K  values, 
the  4>n(K)  is  not  a  decreasing  function,  but  rather  it  in¬ 
creases  to  a  local  maximum  and  then  goes  to  zero  as  the 
modified  Von  Karman  spectrum. 

D.  APPROXIMATE  FORMS  OF  THE  SCALAR  EQUATION 

The  most  coimnonly  used  approximations  for  the  scalar 
equation  are  the  Rytov  (small  or  smooth  perturbation)  approxi¬ 
mation,  and  the  parabolic  equation  approximation. 

The  Rytov  approximation  is  based  on  the  perturbation 
theory  and  has  one  basic  assumption: 


IVlnEj  «  |  V  in  E0  | 


(11-10) 
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This  restriction  is  equivalent  to  th.  assumption  that  the  been 
cin  be  scattered  by  turbulence  eddies  only  once,  wbicb  »»• 
that  it  restrict,  the  range  of  applicability  of  this  appro— 

mation  to  weak  turbulence  levels. 

Starting  with  the  scalar  wave  equation,  and  replacing  n 

by  l+nr  and  E  by  exp(T)  ,  we  get 

V2T  +  (VT)2  +  ^(l  +  fti)2  *  0 


Taking  T  -  T„  ♦  T,  and  n2  «  2n,,  w.  9«t 

,2,  +  2  +  k2  +  ,2^  *  (,V  2  *  2  ,«„)  (VT,)  ♦  2k2n, 


to 


We  can  see 


that  condition  (H-10>  i*  sqniv^ent 
|VT,I  «  l7Tol 

„d  that  if  eT°  -  *0  “  «“  ™  '0lUti°n  °f  ^  SOal“ 
equation*  the  last  equation  can  be  written  a. 

72T,  *  27T, ' VT,  ♦  2kJn,  ■ 


2.  ,  o  «*'“■> 


which  can  be  solved  (VT,  is  known, ,  for  the  — oi  V 
Experimentally  it  wee  found  ,7,  that  th.  kytov  apprckiaa- 

tion  is  valid  as  long  as 


„2  .  .12*  „V.  tW*  C2  <  .3 


(11-12) 
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where : 


wave  number. 


L  ■  distance  of  propagation,  and 


the  turbulence  strength. 


The  Markov  approximation  is  based  on  the  assumption  that 
(for  an  CM  field  which  propagates  along  the  x  axis) 


(11-13) 


where 


u  =  u(x,y,z)  = 


E  (x,y,z) 
exp(lkx) 


(11-14) 


As  Tatarski  showed  [2],  the  condition  (11-13)  is  equivalent 
to  the  condition  1  <<  ki,Q  which  is  almost  always  met  in  the 
atmosphere  (lQ  -  1mm)  for  the  visible  band. 

Inserting  (11-14)  in  the  scalar  equation  and  using  condi¬ 
tion  (11-13)  gives  the  following  equation,  which  is  called 
the  parabolic  equation 


2ik  +  7r  u(x'£oJ  +  2k2n^ (x,^) u(x,£q) 


0  (11-15) 


where  r^  is  a  vector  in  the  y,z  plane. 

The  physical  meaning  of  condition  (11-13)  is  that  the 
scattering  due  to  the  eddies,  is  primarily  in  the  forward 


direction  [8] .  Eq.  (11-15)  can  be  solved  by  using  the  Markov 
approximation  [2]  or  by  path  integral.  Dashen,  who  was  the 
first  to  show  that  the  parabolic  equation  could  be  solved 
by  the  path  integral  technique,  argues  [10]  that  the  con¬ 
straints  imposed  on  the  solution  of  the  parabolic  equation  by 
both  techniques  are  the  same,  and  identical  to  the  one  imposed 
by  the  parabolic  equation. 

E.  PATH  INTEGRALS  AND  THE  SOLUTION  OF  THE  PARABOLIC  EQUATION 

In  the  forties ,  the  path  integral  technique  was  developed 

R.F.  Feynman  in  an  effort  to  use  the  principle  of  least 
action  as  the  basic  principle  of  Quantum  Mechanics.  A  good 
description  of  the  technique  and  its  application  in  Quantum 
Mechanics  and  a  few  other  fields  of  physics  can  be  found  in  a 
textbook  by  Feynman  and  Hibbs  [21] .  The  use  of  the  technique 
in  solving  the  parabolic  equation  is  described  in  a  report 
by  D.L.  Palmer  [22] . 

The  use  of  the  path  integral  technique  in  solving  differ¬ 
ential  equations  can  be  demonstrated  by  the  following  example 
(free  particle  in  a  one  dimension  motion) . 

Applying  the  Shrodinger  equation  to  this  case  gives  us 

2  2 

-ih  ||r  »  ^  (+  boundary  conditions)  (11-16) 

where  <(i  is  the  probability  amplitude  function. 

Me  know  that  for  a  free  particle  in  a  one  dimension 
motion, 
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L  = 


mx 


~T 


and  the  action  integral  is  given  by 

*b 

S  *  /  L  dt  (11-17) 

fca 

Starting  from  Feynman's  basic  assumption,  namely  that 
the  probability  amplitude  assigned  to  the  ith  path  (connecting 
points  a,b)  is 


b 

ij^Ufb)  *  /  exp(is(a,b)/h]dfi 

a 

where  the  integral  is  being  carried  out  along  the  ith  path. 

Since  the  total  probability  amplitude  is  the  sum  of  the 
contributions  from  all  possible  paths,  it  is  given  by 

\j»(a,b)  *  /  D  (paths)  exp  [is  (a,  b)  /h]  (11-18) 

Compar ing  Eqs.  (11-16)  and  (11-18),  we  can  see  that  the 
solution  of  the  differential  equation  (11-16)  can  be  found 
by  the  use  of  the  path  integral  technique,  providing  we  can 
cone  up  with  the  kernel  of  the  differential  equation,  which 
in  our  example  was  obvious. 

In  the  case  of  the  parabolic  equation,  the  procedure  is 
aore  elaborate,  and  since  this  is  not  the  main  topic  of  our 


21 


work,  we  would  rather  merely  quote  the  result  in  Section 
(I1I.C)  of  this  work. 

F.  SCINTILLATION  AND  SATURATION 

Viewing  a  source  of  constant  output  through  the  atmos¬ 
phere,  we  see  that  the  intensity  varies  with  time  (scintil¬ 
lates)  ,  just  as  in  the  case  of  stars  twinkling. 

Usually  the  scintillation  strength  is  measured  by  the 
variance  of  the  normalized  intensity 

a2  «  <(I/I  )2>  -  <I/IQ>2  -  -1  (11-19) 

X  u  u  <I>* 

where : 


Since  it  was  found  that  in  weak  turbulence,  the  log  ampli¬ 
tude  of  the  EM  field  is  a  random  variable  with  normal  distri¬ 
bution  [2] ,  the  log  of  the  normalized  intensity  is  distributed 
normally  in  this  region.  Hence  it  is  useful  to  work  with 
the  variable  l  ■  ln(I/IQ)  and  to  define  the  turbulence  strength 
by 

a2  »  <Z2>  -  <Z>2  (11-20) 

Defining  &I  *  I  -  IQ,  we  can  see  that  if  AI  <<  IQ ,  since 
<&I>  *  0,  then 
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I  *  ln(I/IQ) 


ln[(AI +  I0)/IQ] 


ln(l  +  ££) 
A0 


AI/I. 


Thus  Eq.  (11-20)  can  be  written  as 


o2  «  <(AI/Iq)2>  -  < AI/Iq>2 


<ai2> 


(11-21) 


There  are  several  models  by  which  the  scintillation  can 
be  explained.  The  simplest  model  which  was  originally  used 
by  Tatarski  [1] ,  describes  the  eddies  in  the  atmosphere  as 
spherical  lenses  whose  index  of  refraction  differs  from  their 
surroundings  by  An  *  .  These  lenses  focus  or  defocus 

(depending  on  the  n1  values)  the  rays  onto  the  target  plane. 

It  was  believed  that  this  picture  gives  good  approxima¬ 
tion  only  in  the  weak  turbulence  region  (not  more  than  one 
scattering  event) ,  until  Yura  [23]  showed  that  it  can  be 
applied  to  the  strong  turbulence  region  also,  if  we  incor¬ 
porate  in  this  picture  the  notion  of  the  transverse  coherence 
length  (rQ)  and  its  effect  on  the  focusing  ability  of  the 
eddies.  Due  to  focusing  and  defocusing  of  the  beam  while 
propagating  through  the  atmosphere,  the  transverse  coherence 
length  of  it  is  gradually  reduced  and  therefore,  the  ability 
of  the  eddies  to  focus  the  beam  further,  is  reduced  too. 

Hence  the  growth  of  the  scintillation  strength  is  saturated 
and  it  should  reach  a  maximum. 
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Other  models  which  are  used,  either  in  detailed  calcula¬ 
tions  or  in  qualitative  explanations,  are  the  phase  screen 
model  [6]  and  the  scattering  model  [24]  .  The  first  describes 
the  effect  of  the  atmosphere  on  the  beam  as  a  distortion  of 
the  phase  front  due  to  different  optical  paths  in  the  turbu¬ 
lence.  In  turn,  the  distorted  phase  front  will  produce  an 
interference  pattern  in  the  target  plane.  The  second  views 
the  EM  fields  at  the  target  plane  as  the  sum  of  fields  which 
were  scattered  at  small  angles  by  off-axis  eddies  and  the  direct 
(unscattered)  field. 

Since  in  all  models,  the  scintillation  strength  depends 
on  the  turbulence  strength  (density  and  effectiveness  of  the 
eddies)  and  the  total  number  of  focusing  or  scattering  events 
along  the  propagation  path,  it  is  found  very  frequently  in 
the  literature  that  the  integrated  turbulence  strength  (the 
effective  turbulence  strength)  is  given  by  the  Rytov  param¬ 
eter  which  was  defined  in  Eq.  (11-12) .  Sometimes  it  is 
denoted  by  8  or  a T. 

It  was  found  experimentally  [7]  that  in  the  weak  turbu¬ 
lence  region  (also  called  the  Rytov  region) ,  in  which  the 

2 

intensity  distribution  is  log  normal,  ox  goes  up  to  0.3.  As 
2 

ax  grows  beyond  this  value,  the  intensity  distribution  is 

close  to  the  K  distribution  [25],  the  scintillation  strength 

2 

reaches  a  maximum  [26]  and  starts  to  decrease.  As  ox  grows 
further,  the  intensity  distribution  is  best  described  by  the 
exponential  distribution  [10,27]  and  normalized  intensity 
variance  approaches  1. 


III.  STATISTICAL  DEPENDENCE  OF  THE  TWO 
PARTS  OF  THE  FOLDED  PATH 

A.  PREVIOUS  WORKS 

As  we  mentioned  in  the  first  chapter  of  this  work,  the 
question  of  whether  the  two  parts  of  the  folded  path  are 
statistically  independent,  is  a  basic  one  when  we  try  to 
describe  the  statistical  behaviour  of  the  EM  field  which 
propagates  along  it.  Once  we  can  assume  statistical  indepen¬ 
dence,  we  can  treat  the  two  parts  separately,  applying  the 
well  known  results  for  the  one-way  propagation  to  each  of 
them.  Since  our  main  object  is  the  scintillation  along  a 
folded  path,  with  a  special  emphasis  on  short  distances  be¬ 
tween  the  source  and  the  detector,  the  first  subject  of  our 
analysis  is  to  try  to  define  a  criterion  for  the  statistical 
dependence  between  the  two  parts. 

Among  the  works  that  we  were  able  to  find  in  the  litera¬ 
ture,  which  analyze  the  folded  path  question,  same  like  the 
theoretical  work  by  Yura  [28]  and  the  experimental  work  by 
Ben-Simon  et  al.,  [29],  assume  statistical  independence  of 
the  two  parts.  Yura  gives  the  condition  for  independence  as 


d  >  Ds,  Dd 


where  d  is  the  distance  source-detector  and  0  ,  D ,  axe  the 

s  d 

aperture  diameters  of  the  source  and  the  detector  respectively. 
Since  in  Yura's  paper  the  target  is  diffusive  (volume  of 
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air  molecules)  and  the  reflected  beam  is  totally  incoherent, 
we  cannot  apply  this  result  to  our  case  (specular  target) . 

In  other  works  [30,31,32]  the  question  of  the  indepen¬ 
dence  of  the  two  parts  is  related  to  the  wavefront  tilt  or 
angle  of  arrival  (AOA)  fluctuations  which  depends  on  the  large 
scale  eddies.  Since  we  are  concerned  with  scintillation 
which  is  caused  primarily  by  the  small  scale  eddies,  we  can¬ 
not  use  their  results. 

Before  we  continue  our  analysis,  we  would  like  to  mention 
without  going  into  detailed  explanation,  that  the  usual  cri¬ 
terion  for  independence  between  two  diverging  paths  (their 
end  points  should  be  separated  by  a  few  isoplanatic  patches) , 
cannot  be  applied  in  our  case,  due  to  the  fact  that  in  a  folded 
path  case  certain  parts  of  it  (near  the  target)  might  be  very 
close  and  certainly  correlated  to  some  degree.  As  we  will 
see  later,  this  correlation  affects  the  scintillation  and 
we  cannot  ignore  it. 

B.  THE  EXPRESSION  FOR  THE  PHASE  VARIANCE 

We  start  our  analysis  with  the  effect  of  the  atmospheric 
turbulence  on  the  phase  of  the  propagating  EM  field. 

The  phase  fluctuations  (S)  relative  to  the  vacuum  value, 
due  to  the  non-uniform  medium,  are  given  by  the  integral 
along  the  optical  path. 


i  > 

i  t 


L 

S  *  k  /  n. (x)dx 
0  1 


(III-l) 
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Since  <n^(x)>  58  0 


<S>  *  <*  /  ni(x)dx>  *  o 


k  /  <n1(x)>dx  =  0 


The  variance  of  the  phase  fluctuation  is 


,2  *  <s2>  -  <S>2 


<s2> 


•2>  =  <  (X  f  n^x)**)**  /  n1(x')dx,)> 

A  ” 


k2  fL  dx  /  dx,<n1(x)n1(x')> 

0  0 


<S  > 


2>  ,  k 2  J  <*  /  dx’Bn(x  x‘) 


(III-2) 


Since  in  the  Kolmogorov 


,  Ix-xM  and  it  is  the  same 


model/  Bn(x,xM  -  Bn(y)  where 
ame  function  regardless  of  the 


origin  of  y 


of  y,  <S2>  can  b9  approximated  by 


<S2>  -  *2L  /  Vy)dy 

-L/2 


provided  L  >>  y  where  y  is  the  range  at  which 
Bn(y)  «  Bn(0). 

As  we  can  see  from  Eq.  (III-2) ,  the  variance  of  the  phase 
is  related  directly  to  the  field  of  the  index  of  refraction. 
This  is  the  reason  we  chose  to  utilize  the  phase  variance  in 
the  criterion  for  statistical  dependence/ independence  between 
the  two  parts  of  the  folded  path. 

In  this  part  of  the  work  we  chose  to  work  with  an  infinite 
plane  mirror  as  a  folding  target  for  two  reasons: 

1.  It  is  easy  to  compare  the  case  of  a  folded  path  (with 
a  plane  mirror)  when  the  angle  between  the  two  parts 
of  the  path  is  29,  to  the  case  of  a  straight  path 
(29  »  tt)  (see  Fig.  1)  .  There  is  no  discontinuity  of 
the  rays  and  there  are  no  singular  points  at  the  target, 
as  in  the  case  of  the  comer  cube. 

2.  When  an  infinite,  or  at  least  very  large  mirror  is  used, 
we  do  not  have  to  take  into  account  the  possibility 
that  part  of  the  beam  will  not  be  reflected.  Thus  at 
this  stage,  we  can  neglect  the  effect  of  beam  wander. 

C.  THE  RELATIONSHIPS  BETWEEN  THE  MEANS  OF  THE  EM  FIELD  IN 

VACUUM,  THE  EM  FIELD  IN  A  NON-UNIFORM  MEDIUM,  AND  VALUES 

OF  THE  PHASE  VARIANCES 

Following  Dashen  [10] ,  we  can  write  an  expression  for  e 
(where  E  *  Re(e)  at  a  point  (L,r)  due  to  a  point  source  at 
(0,^).  This  analytical  signal  is  the  solution  of  the  para¬ 
bolic  equation  [22]  and  it  is  given  by 


£(0,ro,L,r)  *  /d (paths) 

L  1  5 

x  exp{ik  /  [i(r'(x))  -n. (r+evx)]dxl  (III-4) 
0  *  -  xx 


where : 

dr(x) 

£’<*>  -  sr- 

e  *  unit  vector  along  x  axis, 

>*v 

From  Eq.  (III-4)  we  can  get  the  expression  for  the  analytical 
signal  of  the  vacuum  field  by  setting  n1  *  0.  Thus  we  get 


e°(0,r0,L,r)  *  /D(paths) 


x  exp{^£  /  [r’CxJl^dx} 


(III-5) 


Taking  the  ensemble  average  of  Eq.  (III-4)  we  get 


<e> 


Jd  (paths)  > 


(III-6) 


iv  L  ,  L 

x  exp{=£  /  Ir' (x) ]^dx  -  ik  /  n. [r(x)+e  x]dx) 
*  0  o  x  ~  ~x 


/D(paths)  x  expi^p  /  lr'(x)]2dx} 


x  <exp{ik  /  n. (r (x) +e„x) dx}> 
0  A  “ 
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Using  Eq.  (III-l)  we  cam  write 

L 

<exp{ik  /  n. (r(x)+e  x)dx}>  =  <exp(S?)> 

0  x  x  i 

P  P  P 

where  =  is  and  S  is  the  phase  fluctuation  along  the  path 

P.  Since  it  is  generally  accepted  [33,34]  that  the  phase 

fluctuation  in  the  atmosphere  is  a  Gaussian  random  variable 

P  . 

with  zero  mean,  is  also  Gaussiam  with  zero  mean,  and 
therefore  it  satisfies  the  following  formula 

S?  <(S*)2> 

<e  >  *  exp[ - j - ] 

Thus 


(*) 

For  Gaussian  ramdom  variable  a,  with  zero  meam  amd  with 
stamdard  deviation  a 


<ea> 


/  eaf (a, 0) da 


"  /  a  ]  dot 

/2rTa  -*>  2a^ 


-L-  Ida  -  S*El2i^I 

/2iTa  -«•  2a  /3iTa  2a* 


»  exp[o2/2]  /  f(a,o2)da  ■  exp[a2/2] 

where  f(a,x)  is  the  normal  distribution  function  for  meam  »  x. 


30 


r  <(SP)2> 
expt-  — 


1 


Using  Eq.  (HI-2)  we  get 


<exp{ik 


/L  n.  (r(x)+e  x)dx}> 
0  1  '  x 


.  e*p{-£  /L  f  Bnt((r(x)-r(x,)2+U-=c')2 

*  0  o 

Assuming  that  the  optical  paths  in  our  case, 
main  contributors  to  the  path  integral,  satisfy 

(r (x)  -  r(x')  !  <<:  I*  -  *’> 

Eq.  (Ill- 7)  can  be  approximated  by 

<exp(ik  JL  ni(r(x)+exx)dx}> 


«?{-£  /L  /L  Bn(x;x')<taax'> 

*  0  o 


Inserting  Eq*-  (*«-»>.  ««-*>•  in  Eq' 

get 

<S2> 

<c(0,tQ,L,r)>  •  e°(0,r<j,L,r)expt  5“" 1 


(III-7) 


1-/2Idxdx* } 


which  are  the 


(III-8) 


(III-9) 


(III-6) ,  we 


(IIX-10) 
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where  <S£>  is  the  variance  of  the  phase  fluctuations,  due  to 
propagation  along  a  path  of  length  L. 

From  Eq.  (III-10)  we  see  that  the  effect  of  the  turbulent 
atmosphere  on  the  average  of  a  propagating  EM  field,  can  be 
expressed  by  the  variance  of  the  phase  fluctuations  along  the 
geometrical  path. 

It  is  obvious  that  in  vacuum,  apart  from  a  phase  shift 
due  to  the  reflection,  the  fields  E^  and  E2l  as  described  in 
Fig.  1,  are  the  same.  We  can  ask  ourselves  what  will  be  the 
relationship  between  <ep>  and  <e2L>?  Since  in  the  case  of 
<e2L>  ttiere  no  Adding  at  all,  it  is  obvious  that  in  this 
case  there  is  no  statistical  dependence  (relating  to  the 
folding).  Hence,  the  closer  <eF>  to  <e2l>,  the  3-ess  depen¬ 
dence  (statistical  correlation  between  the  two  parts) ,  we 
have. 

Using  Eq.  ( 111-10}  and  the  observation  that  up  to  a  phase 
constant  e®  *  e2L'  we  cone  to  conclusion  that  the  statis¬ 
tical  dependence  between  the  two  parts  of  the  folded  path, 

2  2 

can  be  determined  by  the  relationship  between  <SF>  and  <S2L>; 

the  closer  their  values  are,  the  less  the  two  parts  of  the 

folded  path  are  correlated.  Therefore,  in  the  next  steps  we 
2  2 

will  calculate  S2L  and  Sp  (we  omitted  the  angular  brackets 
for  the  sake  of  convenience) ,  compare  their  values,  and  formu¬ 
late  a  criterion  for  statistical  independence. 

Before  we  go  to  the  next  step  of  the  analysis,  we  have 
to  evaluate  the  effect  of  the  constraint  which  was  formulated 
in  Eq.  (III-8) ,  on  the  applicability  of  our  approach. 
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Since  the  point  (x,r(x))  is  a  point  on  the  optical  path 
through  the  turbulent  medium  (which  is  given  by  the  normals 
to  the  wavefronts) ,  the  constraint  (111-8)  is  equivalent  to 


dr  (x) 
I  dx 


<<  1 


which  means  that  the  normals  to  the  wavefronts  are  almost 
parallel  to  the  x  axis,  which  in  turn  means  that  only  small 
angles  of  scattering  are  allowed.  This  constraint  is  the 
same  as  the  one  for  the  utilization  of  the  parabolic  equa¬ 
tion  [8,10]  (Eq.  (11-13)),  and  both  are  met  for  almost  all 
conceivable  conditions. 

D.  THE  VALUE  OF  S \  IN  TEEMS  OF  sj ,  FOR  THE  VON  KARMAN 
SPECTRUM  U  U 

From  Eq.  (111-2)  we  can  write 


j  L  L 

x  /  dx  /  dx '  B  ( |  x  -  x '  j ) 
0  0 


and 


-  ,  2L  2L 

S,T  -  k*  /  dx  /  dx'  B  (|x  -  x1 1 ) 
2L  0  o'  n 


Working  on  both  integrals  and  using  the  transformation 


Z  »  X  -  X ' 
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y 


X  +  X* 


X' 


1 

I 


2  2L.-Z  0 

£-{  /  dz  /  B  ( |  z  | )  dy  +  / 

2  0  z  L. 


2La+z 

dz  /  Bfl  ( 1  z  t )  dy } 
-z 


dz  Bn ( I z I ) (2Li~2z)  + 


0 

/  dz  Bn ( i z { ) (2Li+2z) 

"Li 


where  Li  is  either  L  or  2L. 

Taking  2  +  -z  in  the  second  term,  we  get 

L. 

S?  *  2k2  /  (L.-z)B  (z)dz 

Li  0  • 

9 

inserting  the  explicit  forms  of  Bn(x)  [2]  and  the  Von  Karman 
spectrum  for  <|>n(K) 


Vx) 


4,  /"  *n(K,xt  x2ajc 


$n  (K,x) 


.033  C*(x) 
(K2  + 


Assuming  c2(x)  =*  constant  and  taking  y  *  we  get 
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VL0 

A  /  dz 

0 


(I^/Lg-z) 


dy  y  Sin(zy) 
(y2+l)il/6 


where : 


A 


.033  23  tt  k2 


and  z  =  x/LQ  * 

Using  formula  2.3(11)  of  Ref.  [35]  we  get 


L/Ln 

,  1  0  1/3 

/  dz(Li/L0-z)zx/JK1/3(z) 


where : 

Ax  »  .033  2S/3  ir3/2  k2  C2  L®/3/r  (11/6) 

and  K17-(z)  =  Modified  Bessel  function  of  the  third 

i/J  kind. 

Solving  the  integral  in  Eq.  (III-ll)  (detailed  solution 
is  given  in  Appendix  A) ,  gives  for  (L^/Lg)  >>  1 

,  A1Tr1/2r(5/6) 

s£t  -  --573 — (VV  (III-12) 

which  is  the  same  as  the  expression  which  is  given  by  Tatarski 

2 

[2].  Thus  we  can  see  that  for  large  (L^/Lg) /  is  approxi¬ 
mately  linear  in  and 
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(III-13) 


E.  THE  EXPRESSION  FOR  S^ 

F 

The  geometry  of  the  folded  path  and  the  notations  which 
we  use  in  this  section  are  described  in  Figure  2.  We  assume 
that  the  mirror  is  large  enough  so  that: 

1.  Diffraction  effects  can  be  neglected. 

2.  There  are  no  fluctuations  in  the  amount  of  energy  being 
intercepted  by  it. 

We  also  assume  that  the  field  at  the  point  D  (E_)  can  be 
found  by  the  path  integral  technique,  using  SMD  as  the  geo¬ 
metrical  optics  path.  It  should  be  pointed  out  that  since 
the  calculation  of  the  E?  value  by  the  path  integral  assumes 

dr  (x) 

| — — |  <<  1  (scattering  at  small  angles) 

the  application  of  the  technique  to  this  type  of  path  might 
be  questionable,  for  the  reason  that  it  includes  one  point  (M) 
at  which  the  above  condition  is  certainly  not  satisfied. 

Our  justification  for  the  use  of  the  path  integral  tech¬ 
nique  along  this  type  of  path,  is  based  upon  two  arguments; 

1.  In  the  case  of  fields  propagating  in  vacuum,  we  know 
that  apart  from  a  phase  constant,  EQ  *  E2l’  Therefore, 
we  can  use  the  path  integral  technique  (Eq.  (III-5) ) 
along  the  path  SMD  by  excluding  the  point  M  from  the 
integration. 
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2.  If  the  results  of  our  analysis  in  the  two  extreme 


cases; 

a.  6*0 

b.  9  -*•  tt/4 

match  the  results  which  are  obtained  in  other  ways, 
the  fact  that  at  the  point  M,  is  infinite,  has 

no  effect  on  the  applicability  of  the  path  integral 
technique  to  our  case . 

At  the  end  of  this  section  we  will  discuss  the  second 
argument . 

As  we  saw  in  Section  III.C,  in  order  to  find  we  have 

F 

2 

to  calculate  SF  by 

sj  -  k2  /  /  dxdx'  B  (|x-  x'|> 

F  SMD  n 

t 

where  the  integration  is  along  the  line  SMD.  Thus 


I*  2L 

S2  »  k 2  /  dxdx'  B  { j x  -  x '  | )  +  k2  /  dxdx1  B  (|x-x'|) 

?  0;  n  LJ  n 


5  L  2L  -  2L 

+  x  /  dx  /  dx'  B  (jx-x'j)  +  k  J  dx 
0  L  n  L 


/  dx'  Bn(|x-x'|) 


*  h  *  l2  *  h  *  r4 


where  integration  on  the  interval  0,L  means  integration  along 
SM,  and  on  the  interval  L,2L  means  integration  along  MD. 

In  integrals  I3  and  I4,  x'  is  the  distance  of  the  point 
(yjX^)  from  the  point  M  along  the  line  MD,  and  the  notation 
|x-x'|  stands  for  the  actual  distance  between  the  points, 
which  is  denoted  by  P. 

Using  the  transformation  y  =  x-L?  y'  =  x'-L 

2  L 

I2  =  k  /  /  dydy '  Ba(|y-y'|)  *  Ix 

Since  Ba(|x-x'|)  is  symmetric  in  x,  x',  I  *  1^,  and 
therefore 


0 


L 

/  /  dxdx '  Bn  { |  x  -  x  •  j ) 


(III-14) 


7  L  2L 

+  Jc  j  dx  J  dx'  B  (|x-x'j)  =*  I,  +  I 

0  L  n  -13 


Prom  Eq.  (III-2) 


The  evaluation  of  1^ 

From  Fig.  2  we  can  see  that 
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X' 


(III-16) 


lx  -  X' 


(L-x1)/Cos(20)  ;  dx'  *  -  dXj/Cos(29> 

y  =  (L-x1)Tan(28) 

=  p  =  t(L-x152Tar2(29)  +  (Xj-x)2}172  (IH-17) 


Hence  I3  can  be  written  as 


v2  L 

k  —  f  dx 


■3  -  COS  (26)  2Lsin% 


/  dxx  B^Ptx^)) 


Applying  the  transformation  (x^)  -  (P,o)  and  observing  that 


v  P  Sin(28+a) 
1-  x  -  L - sint  2  0l 


2.  x. 


L  -  P  Sin  a  Cot*29) 


3*  J(P75T) 


3x 

3P 


3 x 
3a 


3xx  3x} 
3P~  la" 


P  cot (20) 


4. 


For  0  <_  a  <_  it/2  -  9 


0  <  P  < 


L  Sin(28) 
Sin(28+a) 


5. 


For  tt/2-9  <_  a  ir-29 


„  ^  L  Sin (28) 

0  I  ?  <  Sin  a 
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we  get 


2  */2-8  *1 

*3  ’  SlW(0^  J  *VP>dP 

tr-26  P2 

+  /  da  J  PB  (P)  dP} 

it/2-9  0 

Making  the  transformation  6  -  ir-(20+a)  in  the  second 
integral,  we  can  see  that  it  is  equal  to  the  first  one.  Thus 

2  tt/2-9  P1 

r3  *  rJtmr  J  **  J  pVp)dp 

where 


p  .  L  Sin(26) 
P1  sinHfl+a) 


Inserting  the  explicit  forms  of  Bn(P)  and  the  Von  Kannan 
$n(K),  and  assuming  that  <j>n(K,x)  »  $n(K),  we  get 


I,  •  A  /  da  /  1  dy  /  dz  z  (z2+l)  ~11/,6sin(yz) 

J  0  0  0 


.033  8ir  k2  C2  lJ/3 
S ini 2© 5 
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•yr— 


where: 


y 


z 


L0K 


p/Lo 


-  L  Sin(26) 
Y1  “  LQSin(29+o) 


Using  formula  2.3(11)  of  Ref.  [35],  we  get 


if/2-e  yi  ... 

I3  *  A.  /  da  /  dy  y  K./3(y) 

0  0  A/'* 


where: 

.033  25/3  ,3/2  k2  C2  I*/3 

A  a  ^  0 

A1  Sin  (29)  fdl/6) 


K  *  Modified  Bessel  function  of  the  third  kind 


f  =  the  gamma  function 


Comparing  the  last  equation  with  Eq.  (III-12)  gives 


22/3  L, 


Sin (29) 


7^ 


( 5/  6 )  L 


tt/2-9 
/  da 


dy  y4/3  K1/3(y) 


Inserting  the  last  equation  and  Eq.  (III-15)  into  Eq. 
(III-14)  gives 


Sp  -  2sJ(l  .  f(y0,e>] 

or  by  using  (III-13) 


43 


(III-18) 


SP  =  S2lH  +  f<V9” 


where 


A-  it/2-8  yl  ... 

f(y0'9)  -  T'  /  da  /  y  '  K1/3(y)dy  (III-19) 

-*0  0  0  J 


and 


A 


2 


y0 


22/3 

ffi/2r (S/6) 
L  Sin (28 ) 


yx  *  y0/Sin(28+a) 


We  solved  Eq.  (III-19)  numerically  and  found  out  that  for 
small  angles  8>  which  is  the  case  in  most  of  the  applications, 
the  f  dependence  on  8  is  very  weak,  and  therefore 
f (Y0»8)  *  f (yQ) . 

The  result  of  the  numerical  integration  is  described  in 
Figure  3,  and  we  can  see  that 


lim  f(y0,8)  »  1  (when  8  0,  yQ  -*■  0) 

y<r° 

lim  f(yn,9)  -  0 

8 -►it/4  0 


and  therefore 


PATHS  DEPENDENCY  FUNCTION 


sJ(L,e) 


i-i 


9  =  tt/4 


(III-20) 


2 

These  are  the  values  of  SF  for  the  two  extreme  points  of 
the  transition  region  (from  exact  folding  to  statistical 
independence  of  the  two  parts  of  the  folded  path) .  In 
previous  works  on  the  subject  of  scintillation  on  folded 
paths,  the  authors  assumed  either  exact  folding  or  statistical 
independence  of  the  two  parts  of  the  folded  paths.  Thus,  as 
we  pointed  out  in  Section  III-E,  we  can  compare  our  results  to 
those  calculated  by  others,  only  at  these  two  points. 

F.  A  CRITERION  FOR  STATISTICAL  INDEPENDENCE 

We  have  now  enough  information  about  the  behaviour  of 

2 

Sp (L, 9 ) ,  in  order  to  define  a  practical  criterion  for  statis¬ 
tical  independence  of  the  two  parts  of  the  folded  path. 
Strictly  speaking,  the  two  parts  are  independent  only  for 
9  *  tt/4.  Yet  because  f(yQ,6)  decreases  very  fast  at  small 
yQ  and  has  a  very  long  tail  (approaches  slowly  to  zero) ,  we 
can  define  the  region  at  which  we  have  statistical  indepen¬ 
dence,  as  the  region  in  which  f(y0,9)  <  e-1.  From  Figure  3 
we  see  that  this  criterion  is  met  (for  small  9) . 

y0  -  4 


Since 


6 


L  Sin (28) 


for  small  angles  8/  yn  can  be  written  as 


yn  -  d/l, 


where  D  is  the  distance  source-detector.  Thus  we  see  that  the 
statistical  independence  between  the  two  parts  of  the  folded 
path  depends  on  the  distance  source-detector,  expressed  in 
units  of  the  outer  scale. 

As  we  said  at  the  beginning  of  this  section,  one  of  the 
justifications  for  employing  the  path  integral  technique  in 
our  case,  even  though  the  main  constraint  on  its  utilization 
is  violated  at  the  point  M,  is  the  observation  that  the  re¬ 
sults  in  the  extreme  cases  (8  =*  0;  9  •  it/4)  match  the  results 
which  are  obtained  by  other  techniques. 

For  the  case  of  8  *  u/4  we  got 

s£(9  -  ,/«>  -  S^L 

which  is  the  expected  result  based  on  physical  arguments  and 
the  results  of  Section  III.D. 

For  the  case  of  8  *  0  we  got 

Sp(e  -  0)  -  2S2t  -  4s£  -  (2St)2 
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When  a  beam  propagates  back  and  forth  along  the  same  turbu¬ 
lent  medium,  the  phase  fluctuation  due  to  each  path  element 
Ax  is  doubled.  Thus  comparing  the  random  variable  S  in  the 
case  of  a  path  of  length  L  (SL)  to  the  case  of  a  folded  path 
with  9*0  (SF) ,  we  get 


=  2S, 


From  probability  theory,  this  gives  for  the  variances 


which  fits  our  results. 


IV.  SCINTILLATION  STRENGTH  FOR  THE  FOLDED  PATH 


A.  THE  RATIO  BETWEEN  THE  SCINTILLATION  STRENGTHS  FOR  THE 
FOLDED  PATH  AND  THE  ONE  WAY  PATH 

The  log  intensity  variance  for  the  Rytov  region  for 
spherical  waves  is  given  by  [1] 


o£2<L)  =  .5  k7'6  C l  L11/6 


(IV-1) 


In  the  strong  turbulence  region  (exponential  distribution 
of  I) ,  the  scintillation  strength  is  given  by  the  variance 
of  the  normalized  intensity  [8] 


«  i  +  -4L 

1  <v' 


(IV-2) 


Comparing  Eqs.  (III-12)  and  (IV-1)  we  get  for  the  Rytov 


region 


6.31 


A(k,LQ)L5/6S2  (IV-3) 


or  for  a  path  of  length  2L 


°l <2l> 


A  a5''6  S221 


(IV-4) 


From  Eg.  (HI-18)  we  get  for  the  folded  path 


r?<F)  -  A  211/6  l5/6  S?(l  +  f(y„,e)l 


We  can  express  the  ratio  between  the  folded  path  scintilla- 
2 

tion  strength  (o^fF))  and  the  one  way  scintillation  strength 
2 

(a  (L) ) ,  in  the  weak  turbulence  region,  by 

co  (F)  11/6 

R  *  -4 -  »  2li/b[l  +  f(yft,6)]  ( IV- 5 ) 

of(L) 

This  ratio  is  the  parameter  which  we  are  going  to  measure  in 
the  experimental  part,  in  order  to  verify  our  model.  The 
dependence  of  R  on  yQ  (for  small  9)  is  described  in  Figure  4. 

B.  DISCUSSION  AND  OBSERVATIONS 

1.  The  Range  of  R  and  the  Saturation  of  Scintillation 
Strength 

As  we  can  see  from  Eq.  (IV- 5)  and  the  corresponding 
figure,  the  range  of  R  is 

211/6  <  R  <  21?/6 

Therefore  when  we  work  with  a  system  for  which  8  is  close  to 
zero  or  zero,  we  should  expect  the  scintillation  strength  to 
be  about  seven  times  greater  than  the  scintillation  strength 
of  the  one-way  path. 

One  should  be  aware  of  this  large  ratio  especially 
while  designing  a  system  or  an  experiment  which  is  based  on 
the  log-normal  distribution.  Since  the  region  of  the  log¬ 
normal  distribution,  as  we  mentioned  in  Chapter  II,  is  given 
by 
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in  the  case  of  exact  folding  (0  =  0) ,  the  upper  limit  of 
this  region  can  be  reached  along  relatively  short  distances 
and/or  relatively  weak  turbulence  levels. 

On  the  other  hand,  since  saturation  can  be  achieved  easily 
(relative  to  the  one  way  case) ,  it  provides  a  way  of  carrying 
out  experiments  in  the  saturated  region  without  the  difficul¬ 
ties  of  using  very  long  paths. 

2.  Comparison  with  Other  Works 

Smith  [36,37]  found  that  the  scintillation  strength 
of  a  spherical  wave  which  propagates  along  a  folded  path 
(plane  mirror)  with  9  *  0  is 

L  oo  2 

a2( F)  =  64  7r2k2  /  dx  /  K$n  (K)  Sin2  [K~* jk£~X)  J  d*  (IV-6) 

Since  he  used  the  phase  screen  approach  [6]  and  the  Kolmogorov 
spectrum,  his  results  are  applicable  only  in  the  Rytov  region 
and  there  is  no  dependence  on  the  outer  scale. 

Using  the  Kolmogorov  spectrum  and  the  Kellin  trans¬ 
form  (formula  6.5(15)  of  Ref.  [35],  we  get  from  Eq.  (IV-6) 

o2  (F,  8  ■  0)  -  A,  /L  dxtx(2fr"*?  ]5/6 
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where 


A,  =  .033(4)  27/6  r(-5/6)(-l)  Cos  (  5tt/12)  it2  k7/6  C2 

j  n 

Making  the  transformation 

y  =  x/L 

we  get 

of{ F,  9  «  0)  *  A_  L11/6  J1  dy  [y  (2-y)  ] 5/6  (IV-7) 

3  0 

It  is  well  known  that  for  the  one  way  path  of  length 

L  [2] 

of  (L)  *  16  ir2k2  /L  d x  /  K*  (K)  Sin2  ]  dK 

x,  0  0  n  * 

Following  the  same  steps  as  in  the  previous  case,  we  get 
o2(L)  »  A4  dx  5/6 

where 

A-  •  .033  4  (-1)  T(-5/6)  Cos  (5tt/12)  ir2  k7/6  C2 

4  n 
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and  we  can  see  that 


Taking  y  -  x/L  we  get 


R 


2  45/6 


1/2  5/6 
I  tx(l-x)]5/6dx 

0 


217/6 


which  is  exactly  the  sane  results  as  we  got  by  using  the  path 
integral  technique. 

It  is  easily  seen  that  even  though  the  results  for 
9-0  are  the  sane  (using  the  two  techniques) ,  Smith's  paper 
cannot  give  the  values  of  R  for  9^0. 
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Game  et  el.,  [38]  measured  the  value  of  R  (for  8  *  0) 

along  a  pat's  of  274m  at  night  (using  a  collimated  beam  HeNe 

laser),  to  be  6.668.  From  their  paper,  it  is  not  clear 

2  2 

whether  the  measurements  of  (L)  and  o ^ (F)  were  done 
simultaneously . 

3.  The  Weight  Function  in  the  Case  of  Exact  Folding  (8  =  0) 

The  relative  contribution  to  the  scintillation  of 
different  points  along  the  path  is  described  by  the  weight 
function. 

From  Eqs.  (rv-7)  and  (IV- 8)  we  can  see  that  the  weight 
functions  VT  (x)  for  spherical  waves  are: 

For  one  way  W_  (x)  *  [x(l-x)]5,/6 

L 

For  the  folded  path  Wp(x)  *  [x(2-x)]5/^6 

The  two  functions  are  shown  in  Figure  5. 

From  Figure  5  we  can  see  that  in  the  case  of  the  folded 

path  (8  =  0) ,  the  main  contribution  to  the  scintillation  comes 

from  the  vicinity  of  the  mirror.  In  the  case  of  the  single 

path, „ the  main  contribution  is  from  the  center  of  the  path. 

Since  it  is  very  unlikely  that  we  cam  find  the  proper  condi- 

2 

tions  under  which  Cn(x)  along  an  atmospheric  path  is  constant, 
we  should  be  careful  when  we  compare  results  from  the  two  types 
of  path. 

< 
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4 .  The  Effect  of  Path  Folding  on  Aperture  Averaging 

The  discussion  and  the  evaluation  of  the  scintillation 
strength  in  previous  sections  assumed  that  we  have  a  point 
detector.  The  term  point  detector  refers  to  a  detector  with 
aperture  diameter  (D^)  smaller  than  a  critical  length  (DQ) . 

The  aperture  averaging  can  be  explained  qualitatively 
in  the  following  way  [7,40].  Due  to  atmospheric  turbulence, 
the  EM  field  at  the  detector  plane  has  a  certain  spatial  struc 
ture.  This  implies  that  the  intensity  (the  modulus  of  the  EM 
field)  has  a  spatial  structure  too.  In  order  to  measure  the 
full  effect  of  the  scintillation  utilizing  direct  detection 
of  intensity,  the  detector  aperture  should  be  smaller  than 
the  distance  along  which  the  intensity  variations  are  negligi¬ 
ble  (Dq) ,  otherwise  the  detector  will  average  different 
intensities  across  its  aperture.  (In  the  language  of  Young 
[40] ,  the  detector  will  filter  out  the  high  spatial  frequen¬ 
cies  of  the  intensity's  spatial  variations.)  The  critical 
length  is  given  by 


BX(D0)  -  Bx(0)  e  A 

where  Bx  is  the  covariance  function  of  the  intensity  at  the 
detector  plane. 

It  is  well  known  [1,2]  that  in  very  weak  turbulence 


D0  'V 


UL) 


1/2 


where  \  is  wavelength  and  L  is  distance  of  propagation. 
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In  the  case  of  strong  turbulence,  it  is  not  clear 
what  is  the  exact  value  of  Oq  [8].  The  difficulty  in  obtain¬ 
ing  the  Dq  value  arises  from  the  difficulties  in  calculating 
the  fourth  moment  of  the  EM  field.  Since  the  transverse 
correlation  length  of  the  EM  field  (rQ)  [23]  is  a  decreasing 
function  of  the  distance  whilst  is  an  increasing  one,  at 
a  certain  distance  the  field  and  the  intensity  spatial  covari¬ 
ance  functions  are  dominated  by  rQ.  Thus,  beyond  certain 
turbulence  strength,  DQ  should  be  a  function  of  rQ. 

As  Yura  [28]  pointed  out,  the  MCF  of  the  EM  field 
is  proportional  to  the  phase  variance. 

It  follows  from  the  above  mentioned  arguments  that, 
as  the  phase  variance  increases  the  critical  length  Dg 
decreases.  Since  we  found  that  in  the  case  of  the  folded 
path 

sj  -  211/6[1  +  f(yQ,0)]  S l 

it  is  obvious  that  while  working  with  the  folded  path,  we 
should  utilize  detector  apertures  which  are  much  smaller 
than  those  which  are  used  in  the  one  way  case. 
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V.  EXPERIMENTAL  WORK 


The  main  object  of  the  experimental  part  is  to  measure 

the  values  of  the  function  R(yQ)  and  to  compare  them  with  the 

analytical  function  which  is  given  in  Eq.  (IV-5) . 

In  order  to  obtain  the  values  of  R(yg),  we  have  to 

measure  simultaneously  the  values  of  the  one-way  scintillation 
2 

strength  (°^)  and  the  scintillation  strength  for  the  folded 
2 

path  (Op) .  These  measurements  should  be  done  at  several  source- 
detector  distances  (assuming  that  Lg  is  the  same  for  all 
measurements) . 

A.  THE  EXPERIMENTAL  SYSTEM 

The  experimental  system  consists  of  two  major  parts: 

1)  The  optical  system; 

2)  Data  reduction  system. 

1.  The  Optical  System 

The  optical  system  is  described  in  Figures  6  and  7. 

The  main  problems  which  we  dealt  with  during  the 
planning  and  construction  phases  of  our  work  were: 

a)  Selecting  proper  sources. 

b)  Making  sure  that  the  beams  which  propagate  along  the 
one-way  path  and  the  folded  path  see  the  same  turbulence. 

c)  Obtaining  high  quality  optical  elements  with  rela¬ 
tively  large  apertures  (within  our  budget  and  time 
span  constraints) . 


59 


I\ 

>  \ 

I  \ 

I  \ 

/  2  6  \ 


Source  2 


Pockels  Cell  Modulator 


Polarizer 

Partially  Transmitting 
Plane  Mirror 


ti 


'A 


Choppe 


Source 


Aperture 


o 


Detector 


Figure  6.  The  optical  system  for  measurements  with  0  ?  0, 
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Figure 


7.  The  optical  system  for  8-0. 
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The  sources  of  our  optical  system  should  meet  certain 
requirements : 

a)  Low  noise.  They  should  have  low  values  of  intensity 
variance  in  order  not  to  obscure  the  intensity  fluc¬ 
tuations  due  to  the  atmospheric  turbulence  along 
relatively  short  paths. 

b)  Operation  in  a  pulsed  mode  or  a  simple  way  of 
modulation.  (This  requirement  evolves  from  the 
way  we  were  planning  to  subtract  the  background 
noise.) 

c)  When  passing  through  the  same  turbulence,  the  respec¬ 
tive  measured  scintillation  strengths  should  be  highly 
correlated. 

d)  Synchronized  operation  either  in  pulsed  or  modulated 
modes  of  operation.  This  requirement  makes  it  possible 

9 

to  compare  pairs  of  pulses. 

Technically  it  is  easy  to  employ  a  GaAs  laser  at  the 
plane  mirror  end  (the  one-way  source)  in  a  pulsed  mode,  and 
a  mechanically  chopped  HeNe  laser  at  the  other  end  (folded 
path  source) .  These  sources  in  these  modes  of  operation, 
were  used  before  by  the  Electro-Optics  group  of  the  Physics 
Department  at  the  NPGS  in  previous  experiments  [41] . 
a .  The  Sources 

It  was  found  (assuming  the  sources'  intensities 
are  normally  distributed)  that  the  normal  intensity  variances 
of  the  sources  were: 
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1)  HeNe  laser  (Spectra  Physics  Model  155)  mechanically 

2  -4 

chopped;  Sj  <  5  x  10 

2  —4 

2)  12  elements  array  of  GaAs?  Sj  <  9  x  io 

2  —4 

3)  8x8  elements  of  GaAs;  Sj  <_  5  *  10  .  This  laser 

matrix  which  is  manufactured  by  Laser  Diode  Labs,  was 
equipped  with  a  Lucite  parallelopiped  rod  mixing  ele¬ 
ment,  in  order  to  enhance  its  output  uniformity.  This 
mixer  was  designed  and  built  at  the  NPGS. 

The  conclusions  from  these  measurements  were: 

1)  In  view  of  the  expected  levels  of  scintillation  strength, 
all  three  lasers  can  be  incorporated  as  sources  in  the 
optical  system. 

2)  The  8x8  GaAs  is  a  better  (more  stable)  source  than  the 

12  elements  array  GaAs.  Therefore  we  decided  to  continue 
the  experiments  with  the  8x8  GaAs  laser  only.  In  the 
following  text,  the  term  "GaAs  laser”  refers  to  this 
laser . 

The  correlation  coefficient  between  the  HeNe  and 
the  GaAs  lasers  were  measured  in  a  configuration  which  is 
described  in  Figure  8.  The  GaAs  in  a  pulsed  mode  was  triggered 
by  the  mechanical  chopper  via  a  pulse  generator.  In  this  way 
we  were  able  to  get  interlaced  series  of  pulses  with  a  con¬ 
trolled  time  separation  between  them.  For  each  laser  we 
measured  the  log-intensity  variance  (in  a  way  which  will  be 
described  in  Chapter  V.B)  and  each  pair  of  variances  (one 
for  each  laser)  was  compared.  For  each  run  we  calculated  the 
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HeNe  laser 


I _ I 


GaAs  pulsed  laser 


Figure  8.  The  configuration  of  the  optical  elements  for 
the  measurement  of  correlation  coefficients 


statistical  correlation  coefficient  between  the  two  series 
of  variances'  values.  One  run  of  this  part  of  the  experiment 
is  described  in  Figure  9.  In  this  particular  run  we  got  the 
highest  value  of  correlation  coefficient  (.973).  In  the  other 
runs  the  correlation  coefficient  values  were  much  smaller. 

As  a  result  of  this  measurement ,  we  decided  to 
replace  the  GaAs  laser  with  a  second  CW  HeNe  laser  (similar 
to  the  first  one) ,  which  is  modulated  by  a  Pockels  cell 
modulator.  (We  tried  to  use  a  second  mechanical  chopper  in¬ 
stead  of  the  Pockel  cell  modulator  but  failed,  because  we  were 
unable  to  synchronize  it  with  the  first  one.)  We  found  that 
the  values  of  the  two  series  of  scintillation  strengths  (log- 
intensity  variances)  of  the  two  HeNe  lasers  were  highly 
correlated.  A  typical  run  with  p  »  .996  is  described  in 
Figure  10. 

The  measurements  of  the  correlation  coefficients 
were  carried  out  in  the  basement  corridor  of  Spanagel  Hall 
at  the  NPGS.  (Its  structure  and  turbulent  features  are  des¬ 
cribed  in  Ref.  41.)  The  path  length  was  145  meters  and  the 
time  separation  between  the  pulses  of  the  two  lasers  was 
250-350  us. 

Based  on  the  results  of  these  experiments,  we 
decided  to  continue  our  experiments  with  the  two  modulated 
CW  HeNe  lasers. 

b.  Frozen  Atmosphere 

As  we  mentioned  in  Chapter  II  of  this  paper,  the 
atmosphere  can  be  considered  frozen  for  time  intervals  of 
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SOURCE  VRRIRNCE 


up  to  1  ms.  Since  the  round  trip  time  of  flight  along  the 
folded  path  (when  the  target  is  up  to  10  km  from  the  source) 
is  less  than  70  ys,  the  atmosphere  can  be  taken  to  be  frozen. 

In  our  experiments  we  employ  two  interlaced 
series  of  pulses  with  a  time  separation  (At)  within  each 
pair  of  pulses.  The  time  separation  in  our  system  is  dic¬ 
tated  by  the  equipment  which  was  available  to  us,  and  is  not 
less  than  250  ys. 

We  tested  the  assumption  of  frozen  atmosphere 
by  measuring  the  correlation  coefficients  of  two  series  of 
log-intensity  variances.  The  two  series  of  pulses  originated 
from  one  source  by  using  two  series  of  slits  on  the  mechani¬ 
cal  chopper  for  the  HeNe  laser,  or  by  the  use  of  the  double 
output  signal  of  the  pulse  generator,  for  the  GaAs  laser. 

For  the  "double  pulse"  HeNe  laser  experiment 
with  t  *  270  ys,  the  correlation  coefficient  was  P*  .998. 
This  experiment  which  was  carried  out  in  "the  corridor"  shows 
that  for  time  separation  of  up  to  270  ys,  the  atmosphere  is 
frozen. 

This  experiment  served  also  as  final  checkup  for 
the  equivalency  of  the  two  channels  of  the  data  reduction 
system  (see  Chapter  V.B) . 

c.  Optical  Elements 

We  originally  planned  on  doing  the  experiments 
with  two  targets;  a  plane  mirror  and  a  corner  cube.  The 
original  configuration  of  the  optical  system  was.  dictated 
by  four  major  requirements: 
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1.  High  quality  elements.  To  maintain  the  phase  relation 
between  different  parts  of  the  beam. 

2.  Large  aperture  elements.  To  intercept  a  large  part 

of  the  beam/  in  order  to  avoid  the  effect  of  different 
transmitted/reflected  intensities  on  the  scintillation 
strength . 

3.  The  laser  at  the  target  end  (the  one-way  transmitter) 
cannot  be  placed  behind  the  target,  because  its  beam 
cannot  be  transmitted  through  the  comer  cube. 

4.  The  beam-splitters  in  use  should  not  cause  an  inter¬ 
ference  pattern. 

The  original  configuration  of  the  optical  system 
is  described  in  Figure  11. 

The  beam-splitters  and  the  plane  mirror  target 
are  made  out  of  large  pieces  of  flat  fused  Quartz,  which  were 
coated  with  Al,  to  achieve  the  required  reflectivities,  at 
the  NPGS  facilities.  They  were  checked  for  flatness  by 
interferometric  techniques.  A  full  description  of  the  testing 
and  coating  processes  is  given  in  Reference  42. 

The  interference  patterns  due  to  the  double  reflec¬ 
tion  from  the  two  surfaces  of  the  beam-splitters,  can  be 
eliminated  by  applying  an  anti-reflection  coating  to  one  of 
the  surfaces.  Since  we  were  not  able  to  use  this  technique 
for  various  reasons,  we  used  a  different  approach. 

The  beam-splitter  at  the  detector  end  is  located 
very  close  to  the  detector  (5-10  cm) .  Hence,  by  using  thick 
plates (taking  apart  the  two  secondary  sources)  we  can  eliminate 
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Figure  11.  The  original  configuration  of  the  optical 
system  for  measurements  with  0  -  0.  a 
is  the  Brewster  angle. 


this  problem  when  the  beam  has  a  small  angle  of  divergence. 

This  beam-splitter  thickness  is  1/4",  it  has  an  elliptical 
shape  (4  x  6"),  and  was  coated  to  give  a  ratio  of  1:2.  The 
beam-splitter  at  the  target  end  is  a  circular  (6”  diameter) 
plate  of  1/2”  thickness.  As  we  can  see  from  Figure  11,  its 
interference  pattern  (at  the  detector)  was  eliminated  by 
applying  the  Al  coating  to  the  back  surface  and  working  with 
a  polarizer  and  angle  of  incidence  equal  to  the  Brewster 
angle,  to  eliminate  the  reflection  from  the  front  surface. 

In  the  course  of  the  experiments,  it  was  found 
by  observing  the  beams  at  the  detector  end,  at  night,  that  the 
beam  which  originates  at  is  very  homogeneous  but  the 
homogeneity  of  the  beam  that  originates  at  S1  suffers  a  great 
degradation  due  to  its  double  pass  through  the  Quartz  beam¬ 
splitter  at  the  target  end.  Since  we  use  a  back  surface 
reflection  for  the  first  beam  (it  passes  this  beam-splitter 
twice  too) ,  the  only  explanation  for  the  different  effects  on 
the  two  beams,  is  that  the  inhomogeneities  in  the  beam¬ 
splitter  material  are  of  larger  length  scale  than  the  diameter 
of  the  beam  near  the  laser  output  aperture. 

Because  of  the  beam  quality  degradation,  we  changed 
the  optical  configuration  from  the  one  which  is  described  in 
Figure  11,  to  the  one  which  is  described  in  Figure  6,  uti¬ 
lizing  the  circular  beam-splitter  which  has  70%  reflectivity 
as  the  target  mirror  (with  the  coated  surface  facing  the 
detector  end) . 
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As  a  result  of  the  above  mentioned  change,  we 
abandoned  our  original  plan  to  measure  the  R  values  with  a 
retroref lector  as  a  target. 

d.  Other  Elements 

1)  The  mechanical  chopper  is  constructed  out  of  circular 
metal  plate  which  is  rotated  at  about  3000  rpm  by  an 
electrical  motor.  There  are  20  slits  on  the  circum¬ 
ference  of  the  metal  plate  with  a  ratio  of  1:4  between 
the  slit's  width  and  the  opaque  regions  between  them. 
The  chopper  modulates  the  CW  HeNe  laser  beam  at  a  fre¬ 
quency  of  about  1  kHz  with  "pulse"  width  (laser  on)  of 
about  200  us.  An  LED  (GaAs)  and  a  detector  which  are 
placed  on  opposite  sides  of  the  chopper  plate,  provide 
a  trigger  signal  at  the  same  frequency  as  the  laser 
modulation  for  the  whole  system.  Due  to  this  arrange¬ 
ment,  we  do  not  have  to  use  a  stabilized  chopper. 

2)  The  Pockels  cell  modulator  is  made  by  Coherent  Asso¬ 
ciates  (Model  3003) . 

3)  The  diameter  of  the  aperture  in  front  of  the  detector 
is  3  mm. 

2.  Data  Reduction  System 

The  data  reduction  system  is  described  in  Figure  12. 
It  is  a  modification  of  a  system  which  was  used  by  the  EO 
group  for  measuring  the  scintillation  of  one  source.  This 
system  and  its  mode  of  operation  are  described  in  detail  in 
Refs.  [41,  42]. 
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Detector 

(Silicon  APD--GE  Model  50HS) 


The  elements  which  were  added  by  us  are: 

1.  The  second  channel  (demodulator  and  log  converter). 

2.  The  electronic  switch. 

3 .  Clock . 

4.  Software  (for  the  Nic-80) . 

The  electronic  switch  (AD  7510D1)  alternately  sends 
the  signals  from  the  parallel  channels  to  the  computer. 

The  clock  is  constructed  of  a  set  of  pulse  generators 
which  are  triggered  by  the  mechanical  chopper's  LED-detector 
pair.  The  time  sequence  of  the  signals  is  described  in 
Figure  13. 

The  change  in  the  software  enables  us  to  store  the 
measured  log-intensity  values  of  the  two  series  in  different 
storing  spaces,  for  later  calculation  of  the  log-intensity 
variances . 

The  main  steps  which  are  carried  out  in  each  channel 
of  the  data  reduction  system,  axe  the  following: 

1.  The  Demodulator  samples  the  signal  intensity  and  later 
(between  pulses)  it  samples  the  background  intensity. 

It  subtracts  the  second  from  the  first  and  sends  this 
difference  to  the  log-converter. 

2.  The  log-converter  takes  the  log  of  the  input  signal  and 
sends  it  to  the  computer  input  via  the  electronic  switch. 

3.  The  computer  samples  the  input  signal,  converts  it  from 
analog  to  digital  and  adds  one  to  the  proper  counter. 
After  sampling  16,384  pulses,  it  calculates  the  variance 
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Figure  13.  Time  sequence  of  the  system's  signals 
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assuming  that  tha  distribution  (aftar  tha  conversion 
to  log  intensity  was  made  by  tha  log  convertor)  is 
normal.  Tha  calculated  value  is  printed  on  the 
assigned  printer. 

Note:  When  we  measured  the  sources'  (lasers')  normal 
intensity  variances,  we  bypassed  the  log  con¬ 
vertors,  and  used  the  average  value  of  the 
intensity  distribution  to  normalize  it. 

After  we  completed  the  modifications,  we  tested  the 
two  parallel  channels  for  equality  by  comparing  two  series  of 
pulses  which  came  from  the  same  laser.  The  high  correlation 
coefficient  between  the  variances  of  the  two  series  (p  *  .998) 
proves  that  the  two  channels  are  equal.  The  results  of  this 
experiment  are  described  in  Figure  14. 

B.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 

The  final  experiment  was  carried  out  along  a  path  of  62 
meters  on  the  lawns  of  the  NPGS.  The  height  of  the  optical 
path  above  the  ground  was  1.1  meters. 

Because  of  weather  conditions,  we  made  measurements  on 
February  27th  and  on  March  3rd  and  4th,  1982.  Throughout  the 
experiment,  the  weather  was  fair  to  clear.  There  were  light 
winds  (0-3  knots)  from  various  directions,  large  portions  of 
tha  path  ware  shaded  by  high  trees,  and  most  of  tha  time  tha 
grass  was  wet.  Tha  meteorological  parameters  which  were  taken 
from  the  NPGS  meteorological  station  (measured  at  the  NPGS) 
with  other  relevant  parameters  are  summarized  in  Table  I. 


TABLE  I 


Information  about  the  conditions  under 
which  the  final  experiment  was  done 


Distance  Source  R  Day  Time  T°F  RH%  Weather 

Detector  (cm) 


0  8.23  3-4-82  1237-1245  56  63  Clear 


20  6.19  2-27-82  1626-1633  58  80  Scattered 

clouds 


40  6.12  3-3-82  1447-1454  56  67  Scattered 

clouds 


70  4.35  2-27-82  1239-1247  61  68  Clear 


105  3.22  2-27-82  1053-1100  58  69  Scattered 

clouds 


During  all  the  runs,  the  wind  was  between  0-3  knots  from 


various  directions 


Figure  15  shows  the  experimental  results,  and  the  function 
R(yQ,9)  which  is  given  in  Eq.  (IV-5) . 

The  R  value  of  each  experimental  point  is  an  average  of 
16  measurements .  Each  measurement  is  the  ratio  of  the  log- 
intensity  variances  of  the  folded  path  (given  distance  source- 
detector)  to  the  one  way. 

The  Yq  values  of  the  experimental  points  were  calculated 
by 


y0,i  “  Di/L0 

where  is  the  distance  source-detector  (for  the  folded  path) , 
and  Lq  was  taken  as  20  cm,  to  given  an  optimal  fit  with  the 
analytical  expression. 

Even  though  we  have  only  five  experimental  points,  it 

9 

seems  that  the  experimental  results  match  the  analytical 
prediction  quite  well  and  confirm  it. 

The  outer  scale  value  (LQ  -  20  cm)  which  gives  an  optimal 
fit  between  the  experimental  results  and  the  analytical  pre¬ 
diction,  seems  to  be  relatively  low,  for  a  path  which  is  1.1 
meters  above  the  ground  in  an  unstable  surface  layer.  For 
a  stable  surface  layer,  this  is  a  reasonable  value  [44] . 
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VI .  SUMMARY 


Optical  scintillation  along  folded  paths  with  different 
angles  between  the  two  parts  of  the  path,  is  a  part  of  the 
noise  in  many  systems  which  measure  the  intensity  of  a  re¬ 
flected  signal  (e.g.,  a  system  which  measures  pollutant 
concentrations  in  the  atmosphere  via  their  absorption) . 

In  previous  works  the  scintillation  strength  along  these 
types  of  paths  was  calculated  by  assuming  either  that  the 
two  parts  of  the  path  are  statistically  independent  or  by 
working  with  an  "exact  folding"  (6  *  0) . 

In  this  work  we  utilized  the  path  integral  technique  and 
the  Von  Karman  spectrum,  to  formulate  a  criterion  for  the 
statistical  independence  of  the  two  parts  (for  spherical 
waves  and  a  plane  mirror) .  This  criterion  was  given  in  terms 
of  the  distance  source-detector  expressed  in  the  outer  scale 
unit ,  namely 


y0  ■  °/Lo  -  4 

For  the  same  case,  we  developed  an  analytical  expression  for 
the  ratio  (R)  of  the  log-intensity  variances  (folded  path 
to  one  way  path)  as  a  function  of  yQ,  in  the  Rytov  region 

R(y0)  »  211/6[1  +  f(y0)3 


where 


f<V 


yQ  =  2L  (L  =  the  distance  source-target) 

y0  -  4 

yQ  *  0  (exact  folding) 


and  showed  that  this  prediction  gives  for  the  end  points 
(yQ  =  O/Yq  =  2L)  the  same  results  which  were  predicted  in 
other  papers  which  dealt  only  with  these  points.  Our  ap¬ 
proach  gives  an  analytical  prediction  for  intermediate  values 
of  yQ,  which  others  do  not. 

In  order  to  verify  our  predictions,  we  designed  and  put 
together  an  experimental  system  by  modifying  an  existing  data- 
reduction  system,  and  constructing  the  optical  system.  This 
system  enabled  us  to  measure  simultaneously  the  scintillation 
strength  along  a  one-way  and  a  folded  path,  which  to  the 
best  of  our  knowledge,  was  never  done  before. 

We  measured  R(yQ)  and  found  that 

1)  the  maximum  and  the  minimum  values  of  R(yQ)  fit  the 
analytical  prediction  well. 

2)  the  slope  of  the  experimental  value  fits  well  the 
predicted  one,  for  an  outer  scale  value  3  20. 

Based  on  these  results,  we  think  that  our  analytical 
expression  predicts  well  the  scintillation  strength  on  folded 
paths . 
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Equations  7.2(12),  7.2(13)  and  7.5(57)  of  Ref.  [43]  give 


the  following  relationships. 


r  (z)  *  l  (z/2)  2®+r/mir  (m+r+1) 

z  m-0 


(A-2) 


Kr(z) 


7T 

Sin  (nr! 


■II_r(z)  - 


Ir(a)l 


(A- 3) 


where  I  (z)  is  the  modified  Bessel  function  of  the  first 
r 

kind. 


H  (z)  -  l  (z/2)2aH'r+1[r  (m+3/2+r)r  (m+3/2)]"1  (A-4) 

m»0 

For  z  -*•  0,  formulae  (A-2)  and  (A-4)  can  be  approximated  by 

Ir(z)  *  (z/2)r/r (r+1)  (A-5) 

Hr(z)  »  (z/2) r+1 [r (3/2) T (3/2+r) ] -1  (A-6) 

Inserting  (A-3) ,  (A-5)  and  (A-6)  in  Eq.  (A-l) ,  and  calcu¬ 
lating  the  lower  limit  gives 


lim  z[K1/>3(z)H_2^3(z)  +  H3y3  (z)K2^3(z)  ] 


^{Siro7TT[I-l/3  (z)  ”*1/3  (z)  ]H-2/3  (z)  2 


IT  i  . 

y  lim 

*  z-*-0 


&tt/3)  tI-2/3(z)“I2/3(z)  1Hl/3  (z }  2 ^ 


n/31/2  lim{ [ (z/2) "1/3/r (2/3) - (z/2) 1/3/r (4/3) ] 
z-*-0 


X  (z/2) 1/3  [r  (3/2)  r  (5/6)  r1 


+  [  (Z/2)  ~2/3/r  d/3)  -  (z/2)  2/3/r  (5/3)  ] 


X  [ (z/2) 4/3 [r (3/2) r (11/6) 3”1! >z 


V31/2  lim{az-b  z5/3  +  cz5/3  -  dz3) 
z-*-0 


where  a,  b,  c  and  d  are  constants. 
Thus  for  z  ®  0 


z [K^^3 (z) H_2^3 (z)  +  H^^3 (z) K2^3 (z ) ] 


(A- 7) 


For  large  z,  formulae  7.13(7),  7.5(52),  7.5(63)  and  7.13(4) 

_2 

of  Ref.  [43]  give,  up  to  0(z  ), 


Kr(z) 


(ir/2z)1/2e~z[l  +  r(3^24-^  (2Z)"1  +  0(z~2)  I 

r(r-  j) 


(A-8) 


(A-9) 


Hr(z) 


ezz”1//2 


r(|+r) 


«(2z)_1  +  0(z“2)  ] 


i  N-l  r  (m  +  i)  0  , 

-  k  l  (-1)“  - 4 - (z/^)*"2®"1  +  Ofz*'21^1)) 

m*0  r  (r  +  j-  m) 


Using  Eqs.  (A-8)  and  (A-9)  in  the  upper  limit  of  Eq.  (A-l) , 


we  get  (z  *  Lj/L0  >>  D 


z[K1/3(z)H_2/3(z)  +  H1/3(z)K2/3(z)] 


(ttz/2)  1/2e""Z{  [1  +  nrKfftZz)"1  +  0(z~2)] 


x  t«*«-.Y2ti  .  M 

/2tT 


“1  +  0(z“2)] 


lr  F(l) 


-  fir tw(z/2) 


~5/3  +  ou*11/3)]] 


+  u  +  Hi7^'(2z)"1  +  0(z~2)1 


[1  -  LO^fl(2z)-1  +  0(z-2)] 


-  ?'rrw  u/2,"2/3  *  ot=_8/3n  1 

*  <fm  -  Wr +  fW  - » <“>’1 


+  0(z”2)}  -  e“z(az1//2  +  (smaller  powers  of  z)  ] 


where  a  is  a  constant. 


i 


which  for  L^/Lq  ■  i  >  10  is  almost  zero  (less  than  10”3) . 
For  the  lower  limit,  by  use  of  Eqs.  (A-2)  and  (A-5) 
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lim  z4/3K.,,(z) 
z-*-0 


lim[24/3/r (-1/3) -az8/3] 


it  2 

Ln  (4rr/ 


Thus  for  L^/Lq  1  10 


Ax  IT  2 


(A-12) 


We  cam  now  write  the  solution  of  Eq.  (IZI-11) 


\  *  ‘ 


(Li/L0)ir1/2r(5/6)  1  r  (13/6)  17,,  ,  , 

1{ - 1273 - C1  +  7lnf7TT  "  Ts]  (VLo) 


+  °<VV~2)  -  nsfif; 


Ax  1T1/2  r(5/6) 


(Li/L0Hl  +  [?£?{!/(?-■■  -  (A-13) 


■(Lj/Lq)"1  +  OC^/Ljj)’2)} 


i 


where 


ir/2-6 


I5^0 


,8) 


/ 


6 a  yf/3K 


4/3(yl} 


and 


yl 


y0/Sin(28+a) 


Taking 


26 +a 


we  get 


I5(y,r-e) 


ir/2 

J 


d6yl/3K4/3(yl)  + 


tt/2+6 

J 

ir/2 


From 


Yi  .  y0/sin  S 


we  get 


dyl 


yn  Cos  b 

- 5 - 

Sin  8 


d6 


4/3<yl> 
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Therefore,  for  the  first  integral  in  1^ 


dB 


Sin2B 


Cos  B 


and  for  the  second  integral 


dB 


,  2, 


Sin  B 


Cos  B 


dYi 


^0 


Since 


Sin2B 


Cos  0 


Sin2B 


u-sin2B)1/2  yi(y? -yo)1/2‘ 


I5(yo'0) 


ya  y1/3KA/,(y) 


v  1/3 

,  “  r  <•  b  y  K4/3(y) 

/  — 2 T172  dy  +  y0  /  — 2 - TT7J 

(yz  -  Yg) 1/2  0  y0 


<y  -  yn> 


where 


L  Sin (20 )/L, 


CB-3-a) 


'a 


yo/Sin(20)  «  L/Lc 


(B-3-b) 


yb  -  y0/S  in  (it/2-9)  -  y0/cos  9  -  (B-3-c) 

Integrating  by  parts  (making  use  of  formula  7.14(3)  of  Ref. 
[43]  gives 
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yi  y4/3KA/,<y> 


y  y(y  -y0) 

o  u 


Since 


+  J  Sec“1(y/y0)y4/3Ki/3^y> 1 

*0 


Sec”1 (1)  =  0> 


Sec“1(ya/y0)  *  ir/2-20f 


Sec”1(yb/yQ)  ■  9' 


and 


y4/\/3(yl 


has  a  finite  value, 


I  (y.,0)  -  - [ (n/2-291 yl/3*4/3 lya’  +  eyb/J*4/3(yb> 


5  VJ0 


*a  yb 


*0 


/  g(y) dy  -  /  gCy)dy 


where 


g(y) 


Sec"1(y/y0)y‘*/  °Ki/3(y)  * 


4/3, 


(B-4) 
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Hence 


A 

f(y0'9)  “  y^tB  (tt/2-0)  +  (Tr/2-0)y^/3K4/3(ya)  (B-5) 

ya  yb 

+  0yb/3K4/3(yb)  -  /  g(y)dy  -  /  g(y)dy 

yo  yo 

This  the  exact  form  of  fty^,©) . 

From  the  definition  of  y  (Eq.  (B-3-b))  we  can  see  that 

a 

in  all  practical  situations ,  whan  the  path  of  propagation  is 
parallel  to  the  ground  but  not  too  high,  we  can  take  y  >>  1. 

cl 

Using  formula  7.13(7)  of  Ref.  [43],  we  get  for  y  >>  1 

yA/\/iW  -  IlhHiA-m)  (2y>~m  - 

-  (ir/2)1/2y5/6a'1[l  +  (55/72)  y"1  +  0  <y"2)  1 

By  taking  only  the  first  two  terms  in  the  brackets,  the  error 
in  the  value  of  y4//3K4y-j  (y)  is  less  than  1%  for  y  >_  10.  Thus 
for  y  >  10  we  can  take 

y4/3K4/3(y)  *  (ir/2)1/2y5/6e“YU  +  t55/72)y"1] 

<  (21T)1/2y5/6e-y 

Comparing  this  expression  with  B,  we  find 
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B/1000  >  y4/3K4/3<Y>  for  y  >  10 

4/3 

Since  in  Eq.  (B-5) ,  B  and  y&'  K^^ly^  have  the  same 
coefficient  (tt/2-0)  ,  we  can  approximate  this  equation  by 

f(y0'e)  *  ~[B(ir/2-0)  +  ey^/3K4/3(yb)  (B— 6 ) 

ya  yb 

/  g(y)dy  -  /  g(y)dy] 

y0  y0 
Numerical  Evaluation 

The  numerical  evaluation  of  Eq.  (B-6)  was  done  in  the 
Computer  Center  of  the  NPGS,  on  an  IBM  3033.  The  program  is 
written  in  FORTRAN  IV  and  is  given  at  the  end  of  this 
appendix . 

The  important  information  about  the  program  is  listed 
below. 

1)  The  function  f(y0#8)  is  evaluated  with  yQ  as  the 
independent  variable.  The  step  by  which  yQ  is  changed 
is  determined  by  a  constant  (parameter)  of  the  program 

2)  The  angle  0  is  calculated  from  the  equation 

8  -  \  Sin”1  (L^q) 

3)  The  parameter  which  is  given  for  each  run,  is  the 
distance  source-mirror  in  units  of  the  outer  scale. 

Lt  -  L/I.0 
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4)  The  constants  A2,  B  are  constants  of  the  program.  Their 
values,  which  were  calculated  from  Egs.  (III-19)  and 
(B-l) ,  are 

A2  =  .7934  B  *  1.1250 

5)  The  values  of  the  modified  Bessel  function  of  the  third 
kind  Kr(y),  are  calculated  in  different  ways  for  the 
two  regions: 

(a)  0.1  £  y  £  4 

(b)  4  <  y 

In  region  (a)  we  used  formula  7.12(21  of  Ref.  [43], 


Kr(y)  *  /  Cosh  (rt)  exp  [- y  fcosh(t)]dt 


which  can  be  written 


Kr(y)  »  /  Cosh (rt) exp [-y  Cosh(t)]dt  +  R(r,tQ) 


Cosh(t)  »  +  6^)  >  jefc 
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Using 


AP-AtZOlfy  NAVAL  postgraduate  SCHOOL  *owtc»ey  ca 

OPTICAL  SCINTILLATION  ON  FOLDED  PATHS. ( U) 
MAH  82  A  ZE  *C V  I 


P/Q 


P'lCLA'jSlf-  110 


NL 


we  get 


R(r,t.)  <  i  /  (ert  +  e*rt)«tp(-  ^f-) dt 

*0 

Taxing 

v  »  efc  x  etQ,  —  ■  dt  and 

x  *  e  ,  t  x 


we  get 


R(r,V  <  i  r  (*r-1+xlr-1) 


for  r  »  4/3;  r-l  *  1/3  end  -(r+1)  *  -7/ 

Therefore 


R(4/3,xfl)  <  \  /  (x+l)e_xy/2dx 


■i  o  -xfty/2  2  -xQy/2 

|[|(x02/y)e  + 


or 


-*0y/2 


R(4/3,xq)  <  y - (x0  +  1  +  2/y) 


R(l/3,xQ)  <  y  J  a  dx  *  e _ 

xQ  y 

Since  in  Eq.  (B-6) 

yb  i  v 

and  we  calculate  f(yQ,e)  only  for  yQ  >.  .1,  setting 
R(r,x0)  <_  10-7  gives 

xQ  <  450  for  r  ■  4/3,  and 

xQ  <  400  for  r  ■  1/3. 

Thus,  we  conclude  that 
6 

K  (y)  »  /  Cosh (rt) exp [-y  Cosh(t)]dt  (B-8) 

r  0 

+  R(r,6) 


where 


R(r,t)  <  10’7  for  r  ■  4/3  or  1/3,  and  y  >  .1. 


We  calculate  Kr(y)  for  .1  <  y  <  4  by  integrating  Eq. 

,-3 


(B-8)  numerically  in  steps  of  At  »  10  .  For  y  >  4, 

we  use  formula  7.13(7)  of  Ref.  [43] 


N-l 


<h)  -  <„/2y)1/V*(  l  U ( 2y)  ~a.O  (y-") )  (B-9) 


We  take  N  »  10,  and  therefore  for  y  >  4,  the  order  of 
the  absolute  error  is  less  than 


(tt/8)  1/2e"44"10  -  io-8 


6)  For  ya  or  yb  which  are  greater  than  30,  we  use  the 
approximation 


v*b 


30 


/  g(y)dy  a  /  g(y)dy 

)  *0 


(8-10) 


For  y^  >  30  we  can  write 


/l  30 

/  g<y)dy  -  /  g(y)dy  +  R(y.) , 

*o  I'd 


where 


R<yi)  “  3Q/  s«c“1(y/y0)y4/3K1/3(y)dy 


<  (ir/2)  /  y4/V/3(y)dy 

30  ■L/J 


4/3  ^ 

<w/2)l-y4/3K4/  (y)] 

30 


For  y  £  30,  (y)  i«  *  decreasing  function. 
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Thus 


R(yx)  i  (*/2) (yx  -  30)304/3K4/3(30) . 

4 

Taking  30  <  <.  10  , 

R(30  <  yL  <  104)  <  (ir/2)  3/2104305/6e“30  -  10“7 

7)  The  numerical  integration  of  / g(y)dy  is  done  in  steps 
of  Ay  *  .01. 

Additional  information  concerning  the  numerical 
integration  can  be  learnt  from  the  program  itself. 

The  Program 
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F( YOt TE TA )  CALCULATION 


SET  CONSTANTS 

.59265 

_ l$0/2. 

8-1.125040 
YO».l 

«-t8i 


?  YO*  »6X  * ‘TETA*  tlOX* • YB*  ,12X,  •  FIRST  TFRH*  ,  4X.  •  YBKtYBP  .  7 


7X 

4 


CALCULATE  YO, VB»  TFT  A 

Si-K4»lTh300 

iM'4 

TET-ARSIN(Z) 

TETA-TET/2. 

YB»Y0/C0SITFTA! 

CALCULATE  FIRST  TERN 

T  ERF= f °  T2-T  ETA l*B 

CALCULATE  SECOND  TERN 

CALCULATE  KSUBI4/3) 

XN-4./3. 

BL*BSL( XN,YBI 

itemslnsr 

CALCULATE  THIRO  ANO  FORTH  TERNS 


c 

c 

c 

c 


XIA*0. 

XtB-O. 

X-I.Z3. 


ii5vsic?y/?SiJ*Soi?vnfY  I 


YO/YOI-O 


Y*30. 

Z*YO/Y 
XSC*AACOS(Z ) 


itHttwa 


XS*X--  . 

X IA»XIA*XS 
Y*YB 

BL-BSMX.Y) 
XS«TETA*<Y* 
XIB®X!B*X5 
XH»I30.-Y0)/DG 
M®IFIX(XM) 

on  102  »  M 

X J*Ft OAT ( Jl 
Y«YO*DG*XJ 


XN) *BL /?• 


*XNI*BL/2. 


102  ^ 


C 


ARCOSIZ  I 

IkrSlW&iU.xN, 

IFfYUcTtYBI  GO  TO 

XIB-XIB+XS 

CONTINUE 

:^i?y^oG 

JI-YB-YO 


*06 

FOR  0G-.01 

IFIY.LF.2T.  I 

TERFO»XlB*Df 
FT-TERF-TEIiSr 


0E«Z1 


TERT-TERFO 


SfftlSI?K%.T|U,YOtT|?F#TERSi,T»S,TEPT,TERFO,FT,FTeTA 


10 


STOP 

END 


FUNCTION  BSLlRtS) 

R*NUt  S*Y »BSL=KNUIY I  MODIFIED  BESSEL  FUNCTION  OF  TMF  THIRD  KIND 
I F( S.GE.4.0 }  GOTO  1 

DU,  TO  WU 

THE1lNT6GRlflONTISNTRUNCATiD,AT  T-6.0  .ERROR  DUE  TO  TRUNCATION  <10-7 
CUT  *-100 • 

DELTA*. 001 

M*6000 

Ml-M+I 

DO*  103  K*I»  Ml 

XKA*FLOAT(K II 

T»XKA*DELTA 

TN«T*R 

tismui 

CC-C-l,l*AC*§ 

TO  o 

XG«|XP(C1) 

GO  TO  28 

28  Jf7x!2o!i.or.k.eo.mii  XG-XG/2. 
bk-mux£ 

03  W9BHfLT8 


GO  TO  8 
CONTtNUI 
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